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MATHEMATICAL MODEL WITH EXPERIMENTAL VERIFICATION FOR 
THE DYNAMIC BEHAVIOR OF A SINGLE -TUBE CONDENSER 
by Andrew A. Schoenberg 
Lewis Research Center 

SUMMARY 

A mathematical model is derived to predict the dynamic behavior of a single-tube 
condenser. The tube is assumed to contain a mobile liquid-vapor interface and to be 
cooled by a gas or liquid stream in crossflow. The model consists of transfer functions 
relating the condensing pressure and the interface position to various disturbance varia- 
bles. The transfer functions are derived from the linearized form of the fundamental 
heat-transfer and flow equations of the process. 

The part of the model relating the condensing pressure and the interface movement 
to vapor flow variations is verified by comparison of predicted with experimental fre- 
quency responses of a gas- cooled mercury condenser. Good agreement is found for the 
amplitude and the phase of the condensing pressure response and the amplitude response 
of the interface movement. The phase lag of the interface movement, however, is con- 
siderably greater than predicted. This disagreement is attributed to the liquid droplet 
transport delay time which was neglected in the theory. 

The two verified transfer functions are in the form of third- over fourth- order poly- 
nomials. The break frequencies of the transfer functions depend on the steady- state oper- 
ating point. This condition is illustrated by changes in the experimental frequency re- 
sponse of the pressure and interface when the operating point is changed. The same 
change is predicted by the model, thus giving further evidence of its validity. 

An additional investigation is made of the variations of the model due to uncertainty 
in evaluating some of the transfer function parameters. The predicted frequency re- 
sponses are found to be sensitive to several of the process properties, such as the heat- 
transfer coefficients and the vapor and liquid pressure drop functions. 

Some of the values of the coefficients of the transfer functions for the cases investi- 
gated are very small and, hence, have little influence on the predicted frequency re- 
sponses. When these values are neglected, a much simpler form of the transfer functions 
is arrived at, which agrees well at the lower frequencies with the original more complex 
model. 



INTRODUCTION 


The condensation of the working fluid in Rankine power cycles has gained increased 
attention due to the development work on space electric power generators such as 
SNAP-8. In these high performance systems where reliability and unattended operation 
are essential, the excursions of the condensing process from the design conditions due to 
various disturbances must be kept within narrow limits. Changes of the condensing pres- 
sure are of special importance because high turbine back pressure reduces efficiency, 
while low pressure may cause pump cavitation and eventual system deterioration. To 
predict these pressure variations and to control them, the dynamic or transient behavior 
of the condensing process must be understood. 

Previous work in heat exchanger dynamics (refs. 1 and 2) has treated the response 
of liquid coolant temperature due to variations in condensing pressure. However, no 
analytic or experimental studies were found where condensing pressure is assumed to be 
the dependent variable, that is, where the disturbances in the system originate in the 
vapor flow or coolant variations. An additional constraint, which has not been considered 
by previous investigators, is the existence of a mobile liquid-vapor interface. In the type 
of system of interest here, the position of the interface in the exchanger tubes determines 
the heat-transfer area of the condensing region. This has the effect of coupling the con- 
densing pressure to the amount of liquid inventory in the exchanger. Thus, there is a 
need to study the dynamics of this type of condenser, and in particular, the behavior of 
the condensing pressure and the vapor liquid interface. 

Such studies, both experimental and analytical, were undertaken at the NASA Lewis 
Research Center. The experimental part of the program dealt specifically with mercury 
vapor condensing inside single tubes. A constant diameter steel tube cooled by gas in 
crossflow was used for the dynamic studies. (This choice of configuration simplified not 
only the instrumentation and measurements, but also eliminated the complication of cool- 
ant dynamics which must be considered in parallel or counterflow liquid cooled conden- 
ser. ) In support of the experiments, as well as to develop an ability to predict time vari- 
able behavior, a theoretical model of the process was derived. The purpose of this re- 
port is to present this derivation and to show how the theory compares with experimental 
results. 

The approach taken in deriving the theoretical model is to approximate the fundamen- 
tal dynamic phenomena rather than to attempt a rigorous solution of the partial differen- 
tial equations governing the condenser. The equations are simplified and linearized for 
small perturbations about a known steady- state operating point and then Laplace trans- 
formed to obtain the system transfer functions. This method has been used with success 
in the analysis of heat exchanger dynamics by various investigators (refs. 1 to 3). The 
derived transfer functions are then analyzed in some detail in terms of the theoretical 
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frequency response they predict. In this analysis an attempt is made to relate the dy- 
namic characteristics to basic physical properties of the condensing process, such as the 
thermal capacitance of the tube walls, the liquid storage, and the pressure drop gradi- 
ents. 

In the last part of the report the predicted frequency responses for the single-tube 
condenser are compared with the experimentally obtained responses to establish the accu- 
racy of the theory. The experimental results used in the comparison are part of the data 
obtained for the mercury condensing dynamics program. 

Two questions arising in the determination of the theoretical frequency response are 
examined in separate sections of the report. One section deals with the error in predic- 
tion of frequency response due to uncertainty of some of the process properties. The 
other is concerned with the possible simplification of the transfer functions. 

Although the derivation and discussions are applied specifically to a mercury con- 
densing tube, the theory is general enough to be applicable to other condensers which are 
uniformly cooled and contain a vapor-liquid interface. For added generality, the deriva- 
tion of the basic condenser element is carried out separately from the inlet and outlet 
pressure-to-flow characteristic of the condenser. This permits simple modification of 
the model so as to include other condenser inlet and outlet conditions. 

THEORY DERIVATION 

Herein is presented the derivation of the mathematical model used to predict the 
variations of the condensing pressure and interface due to various disturbances in the 
process. The disturbance variables are assumed to be the vapor input flow rate, the re- 
ceiver pressure, and the coolant temperature. 

The analysis itself is divided into three sections. The first deals with the derivation 
of the transfer functions for the condensing pressure variations, assuming the interface 
movement to be an independent variable. In the second section, the transfer functions 
for the interface are derived, assuming in turn that the condensing pressure is an input 
variable. In the last part of the derivation, the two transfer functions are analyzed and 
then combined to obtain the desired model for the pressure and interface behavior of the 
total condensing process. 


Condenser Description 

The condenser that is analyzed is shown schematically in figure 1. It consists of a 
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single tube with inside and outside diameters 
Dj and Dg, respectively. The tube wall 
temperature T,_(x, t) is a function of both the 
axial position from the inlet x and time t. 
The tube is convectively cooled by a cross- 
flow gas stream at a temperature Tj^(t). 

The inlet vapor flow rate Wjj is determined 
by a choked valve so that this flow rate into 
the tube can be assumed independent of any 
changes in the downstream conditions. 

The total condensing flow rate w depends on the heat flow into the wall of the tube. 
This heat flow is proportional to the difference between the condensing temperature 
T (x, t) and the wall temperature T (x, t) multiplied by the coefficient of heat transfer. 
The vapor has to condense completely over the condensing length S.^, which is deter- 
mined by the location of the liquid- vapor interface. The remaining length of the tube 
is filled with liquid that is subcooled and then exhausted into the receiver. 

The liquid outflow rate Wj^ is determined by the liquid flow characteristics of the 
valve at the exit of the tube and the difference in pressures at the interface Pj and the 
receiver p^. 



Figure 1. - Schematic of condensing tube. 


Derivation of Condensing Pressure Transfer Function with 
Interface Position as Independent Variable 


The equations that govern the heat transfer and fluid flow of the condenser just de- 
scribed are 

(1 - Xq)wjj(t) + w^(t) - = -Pl — (1) 

W (2) 

r n 


dT^(x, t) 
dt 




0 c 

w 


(°i - 


[T^(x,t)-T^(x,t)] 


®2»nw 


Pw'w 






[T^(x,t)-Tj^(t)] 


( 4 ) 
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(All symbols are defined in appendix A. ) Equation (1) expresses the continuity relation 
for the liquid phase. It states that the sum of the condensing flow w and the liquid in- 
put flow rate (1 - X^^)Wjj minus the liquid outflow Wj^ must equal the net rate of decrease 
of the product of gas phase volume and the liquid density p-^. The quality of the input 
stream is X^. The second equation similarly states the conservation of mass in the gas 
phase. The density of the gas is for simplicity taken to be that at the tube inlet. 

The third equation represents the heat-transfer process in the vapor region of the 
condenser. The right-hand side is the heat flux into the wall. The left-hand side is the 
heat flux given up by the vapor during condensation. The assumptions made here are that 
(1) the amount of internal energy change in the gas phase is negligible and (2) the heat 
content of the input stream is the heat of vaporization with a negligible amount of 
superheat energy. 

The fourth equation states the heat flux balance from the condensing region to the 
coolant. It includes a heat storage term due to the wall heat capacity. The wall is con- 
sidered as having a uniform temperature in the radial direction. 

Two additional assumptions are made to allow the integration of equation (3). The 

coefficient of heat transfer between condensate and wall h is assumed constant. Its 

w 

variation will generally not affect the heat transfer significantly since its value is much 
larger than that of the coolant side. This is particularly true for the gas-cooled tube 
used in this investigation. 

The other assumption is that the condensing temperature variation along the length 
X can be approximated in the steady state by 

T^.(x, t) = T^(t)(l - mx) (5) 

This form is introduced to account for a significant pressure drop in the vapor region 
noted for some operating conditions of the condenser tube. An example of the tempera- 
ture profile is given in figure 19 in appendix B. The effect of equation (5) is to introduce 
the average of inlet and outlet condensing temperature for the steady-state heat-transfer 
calculations. Variations in m are neglected so that dynamically the average tempera- 
ture is assumed to vary as the inlet temperature. 

To relate the gas phase density and the condensing temperature to the condensing 
pressure, the perfect gas law and the Clausius-Clapeyron equation, respectively, are 
used: 


Pv 



( 6 ) 


and where R = gas constant 
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^Pn H,p, 


L^'c 


dT, 


( 7 ) 


rt: 


Another relation needed for the simultaneous solution of the first four equations is 


Vv<‘)=Vc<‘> 


( 8 ) 


where V (t) is the volume of the gas phase, A the cross-sectional area of the tube, and 

V c 

£ (t) the condensing length. 

Several assumptions are implied in this formulation of the vapor volume. First, it 
is assumed that the liquid droplets inside the vapor take up an insignificant amount of vol- 
ume. This is justified due to the large density difference between the gas and liquid 
phases (3000:1 for mercury at 16 psia). Also assumed is a distinct liquid-vapor interface 
that moves inside the tube in a piston- like manner when the liquid inventory of the conden- 
ser changes. A less obvious implication of equation (8) together with equation (1) is that 
the liquid condensate collects at the interface without significant time delay. This as- 
sumption, although not quite accurate, avoids the great complexity of formulating a valid 
model of droplet formation and velocity. 

Conceptually these eight equations plus specified initial conditions are sufficient to 
determine the condensing pressure in terms of the input flow Wjj, the condensing length, 
and the coolant gas temperature. The nonlinearity of some of the equations and the com- 
plex interdependencies of the variables, however, make it more practical to solve them 
in a simpler linearized form. 

The simplification consists of linearizing the equations by considering only small 
perturbations of the variables about a known steady-state operating point. The linearized 
differential equations can then be Laplace transformed and solved algebraically to obtain 
the transfer functions relating the condensing pressure to the disturbance variables. 

The linearization does limit the predictive ability of the model, particularly when 
large changes in the operating point are considered. Where the disturbances are reason- 
ably small, however, such as expected during normal operation of the Rankine loop or 
for sinusoidal testing, the simplifications will not detract significantly from the accuracy 
of the predicted dynamic response. 

The variation of the variables about their steady- state values can be written as 

Wjj(t) = Wjj + AWjj(t) 
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and similarly for all other variables. Substitutii^ these into equations (1) to (8), elimin- 
ating the steady-state parts, and Laplace transforming with respect to the time variable 
result in the following equations (s is the Laplace operator) : 


(1 - Xq)AWjj(s) + AWj.(s) - AWl(s) = -Pj^s AV^(s) 
Xq AWjj(s) - Aw^.(s) = p^s AV^(s) + V^s Ap^(s) 

Hl Aw^(s) = h^^,rDi[AT^(x, s) - AT^(x, s)]dx 

s)(t^s + 1) = AT^.(x, s) + AT^(s) 


where 




4(hcw^^l + *^nw^2) 




^w^l 


^w^l+\w^2 


^2w- 


h Do 
nw 2 


‘^cw^l + *^nw^2 


AT^,(x,s) = ATp(s)(l - mx) 


Ap (s) AT (s)p 
Ap»=— 


RT, 


RT^ 


dPc 
c (J.J. 


AT^(s) 


AVy(s) = A£^(s) 


(la) 

(2a) 


(3a) 

(4a) 


(5a) 

(6a) 

(7a) 

(8a) 
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The second term on the right-hand side of 
equation (3a) represents the variation of heat flux 
due to changes in condensing length. It is the 
product of the integrand evaluated at the interface 
position i and the increment of interface 

The variations represented by equations (3a) 
and (7a) are shown schematically in figure 2, 

Note the change in overall level of condensing 
temperature AT (t) at a constant slope of the 
profile and the change in condensing length, which is accounted for by the change in the 
limit of integration. The terms of the other equations are obtained directly by proper 
substitutions of the variations. 

The integral of equation (3a) can be evaluated by assuming a constant h^^ and by 
employing the relations of equations (4a) and (5a). In dimensionless form, equation (3a) 
becomes 



Figure 2. - Idealizedtemperature profile for condensing 
region showing expected variations. 


where 


= z, 

C L 


(^Iw® AT^(s) C 


n 




+ L 


W„ = 


Aw^(s) 




mf 


Zl-1- 


T = T Z, - T 
a c L n 




AT (s) 
C = — 

n m 


(9) 
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A£.(s) 


and where the following steady- state equalities have been used: 


w H, = 
c 1 j 


^cw^^l^nw^2^c'^a 
^cw^l ‘^nw^2 


w 


:«L = '>ew’''>l'c[Tc(‘J-Tw(«c)] 


Equation (9) represents the overall dynamic heat balance of the condensing region, which 
is the basis of the condenser model developed in this report. 

The condensing flow W and the condensing temperature variation AT /T can be 
written in terms of the input flow, the condensing length, and the pressure by employing 
equations (2a), (6a), (7a), and (8a). The desired relations for condensing flow is 

= Wg - - (1 - /3^)9SP^ (10) 


where 


W„ = 


AWjj(s) 


H 


w 


H 


= Py'^c^c 

Vh 






in which for steady state 


H " 
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The relation for the condensing temperature variation is 


where again 



( 11 ) 


Substituting these equations into equation (9), combining like terms, and solving for 
yield 


P- = 


K, 




— (1 +0s)L 
T S + 1 ^ ' 


W 


( 12 ) 


where 


1 + 

Gg(s) = 

asS^ + bgS + 1 


T (1 - (3 )6 

= 

K, 


>=s ^ ^Iw + 


(1 - 

K. 


(13) 


and 



dT 


c 


^Pc 


Equation (12) represents the desired dynamic model, which relates the condensing pres- 
sure variations to the inlet vapor flow, the coolant temperature, and the interface varia- 
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Figure 3. - Block diagram for condensing pressure varia- 
tion with condensing length as independent input. 


tions. Figure 3 shows in block diagram form the 
relations involved. 

The model is expressed in the form of transfer 
functions which may be used directly to obtain the 
frequency response characteristic for sinusoidal 
input variables or may be transformed back to the 
time domain for arbitrary disturbance functions of 
time. The variables are expressed in dimension- 
less form to emphasize the basic dynamic elements 
and the essential parameters of the model. The model is sufficient to predict dynamic be- 
havior where the input variables, particularly the interface movement, are specified in- 
dependently. In the case of the condenser tube used for the dynamic tests or condensers 
of similar geometries, however, where the interface movement is not independent of the 
other process variables, the dynamic model is not complete. The interactions between 
the condensing pressure and the interface movement must be established. This is done 
in the next section. 


Derivation of Interface Transfer Function 

The interface movement is strongly dependent on the exit condition and the pressure 
drop characteristics of the condenser. The derivation of the mathematical model de- 
scribing this movement is carried out for the case where the condensate exhausts into a 
receiver, whose pressure is independent of the process. This is not the most general 
exit condition; however, the equations can easily be adapted to fit other cases such as a 
pump. 

The receiver has the special characteristic of controlling the pressure at the inter- 
face. This is a result of the small pressure drop that exists in the liquid column between 
the interface and the receiver. The pressure is controlled by the movement of interface 
due to the imbalance between liquid outflow and condensing flow, which occurs in the di- 
rection so as to restore the pressure equilibrium on the liquid column. 

The equations relating the interface movement to the various flows are already for- 
mulated in the previous section and are expressed in linearized form by equations (la) 
and (2a). Adding the two equations eliminates the unknown condensing flow and results in 
the expression for interface movement as given in equation (14), where the variables are 
again expressed in dimensionless form: 
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+ 


(14) 


L 


c 


Wl-^h 




(1 - ^x)Pc 

ft., 

Py 


Equation (14) can be simplified by noting that the term (1 - /3jj)P^/[(Pj^/p^) - l], which 
accounts for the change in mass in the gas phase due to changes in pressure, is very 
small in comparison to mass changes due to liquid changes in the interface. With the 
pressure term eliminated, equation (14) takes the form 


L 


c 



(15) 


The next step in the derivation is to determine the liquid outflow Wj^ in terms of the 
known process variables. The approach to the problem is to consider the behavior of a 
liquid column moving inside a tube with frictional pressure drop. The driving force for 
this movement is the fluctuation of the pressure difference between the interface and the 
receiver. These pressures are coupled to the condensing inlet pressure through the two- 
phase and the liquid pressure drop characteristic of the condenser tube. 

The fluid force balance for the liquid column inside the tube is 


(Pf - Pr)Ac = A^.qL(wL) + 



(16) 


where (p^ - p^ is the pressure difference between the interface and the receiver, qjj(wj^) 
represents the frictional pressure drop in the liquid region which is assumed a function 
of flow rate only, and Mj^ is the inertial mass of the liquid column being accelerated 
relative to the tube at the rate dvj^/dt. Also the following equalities should be noted: 



(17) 


(18) 

^L^c 



The pressure at the interface, as given by equation (19), differs from the condensing in- 
let pressure by the pressure drop in the vapor region q^. This pressure drop depends 
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primarily on the two-phase friction factor, the velocity, and the length of the vapor re- 
gion. Thus, in general, is assumed to vary with flow rate and condensing length: 

Pf = Pc-qv(wH,^c> 

Substituting equations (17), (18), and (19) into equation (16) and linearizing yield 


9Wh ^ a^c 


dt 


( 20 ) 


where Sq^/9wjj, 9q^/9Lp, and are partials of the pressure drop which must 

be evaluated from steady- state experimental or theoretical relations. 

It should be noted that to be strictly correct, the slope m of the temperature pro- 
file of equation (5) should be coupled to the pressure drop. For simplicity, however, the 
effect of the pressure drop variations on m is neglected. 

Laplace transforming and changing to dimensionless form allow equation (20) to be 
solved for Wj^: 


^ 3<Po - PrLVH_lV:c 

^ TpS + 1 

where 


( 21 ) 


K_ 



P 


r 


APr 


Kw 
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^ diL 


T 


P 



(22a) 


The conversion factor g must be included in equation (2 la) for English units where 
is in psi; therefore, 




A cr 

A„ g 


9w, 


(22b) 


Equation (21) is now combined with equation (15) to eliminate Wj^: 


L„ = 




yPc - - Kw^H - ^r^c 


W, 


1+TpS 


H 


(23) 


where 


^r ^ V>— ■ 


Solving for and combining similar terms give the desired transfer functions relating 
the interface movement to the other condenser variables: 


Lc = G.(s)(P^ - Pp - (^V ^ 


(24) 


where 
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(25) 


G.(s) 


K 


K^(as‘‘ + bs + 1) 


e T 

= rp 





^Ip" 





1 + K„ 




:-k> 

P. 






K^,(l + bs + as^) 



Figure 4. - Block diagram for dynamic model of 
Interface movement. 


Equation (24) is shown in block diagram form in figure 4. 

This completes the derivation of the mathematical 
model for the dynamic behavior of the interface for the 
case where the condensate exhausts into a pressurized 
receiver. In the next section this behavior and how it in- 
teracts with the condensing pressure will be examined. 


Analysis and Synthesis of Transfer Functions 


In this section a closer look will be taken at the two parts of the dynamic model de- 
rived in the previous sections in order to arrive at a better understanding of the elements 
that determine the time- variable behavior of the condensing process. This is done by 
analyzing the frequency response characteristics and the coefficients of the major trans- 
fer functions. In the last part of this section the derivation of the model of the condenser 
as a whole interacting process is completed. This is done by the synthesis of the two 
previously derived elements. 

The major dynamic element of the model shown in figure 3 is the transfer function 

defined by equation (13), G (s). It relates the condensing pressure to the input flow vari- 

s 

ations when the other input variables of equation (12), C and L , are assumed zero. 

n o 

This transfer function in factored form is 





1 + ^w® 


(TjS + 1)(t2S + 1) 


(26) 
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where 



This factoring is convenient when the parameter inside the square root is less 

than 1. This also means that the denominator has two positive real roots or, in control 
terms, the system is overdamped. 

Examination of the constants shows that this overdamped situation prevails as long 
as the time constants of the wall t and are considerably larger than the vapor 
capacity time constant 9 . These time constants have the following physical significance. 
The constant is the time, in seconds, required to raise the average wall tempera- 

ture by 1° when the heat flow into the wall is induced by a constant 1° temperature differ- 
ence between the wall and the coolant. The definition of differs from that of 
only in that the heat flux into the wall is induced by a unit temperature difference between 
the wall and both the coolant and the condensate side. Thus, t is always greater than 
T^. The time constant 9 (see eq. (10)) is the time in seconds for the mass in the gas 
phase PyVy to be replaced at the steady- state vapor input flow rate X^Wg. 

For the case where the parameter ^a^y^g is small, say less than 0. 1, a conceptu- 
ally simplified transfer function may be obtained by expanding the square root term in a 
Taylor series where only the first term is significant: 



Skibstituting this approximation into equations (27) and (28) results in a simplified form of 
Tj and Tgt 


i 
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Phase lag, <p, deg Log magnitude 


(29) 


T 



la 




b 


s 


and 




Iw 


(30) 


The approximate frequency response is obtained by substituting the values of 
and T 2 ^ and the sinusoidal frequency 2jrjf for the operator s in equation (26) (note that 
f is in units of cps) : 


Gg(2jrjf) = 


1 + jfr^ 

(1 + 27TjfTj^)(l + ZffjfTg^) 


(31) 


Since the general inequality relation 


^2a>^w>'^la 



is known, the shape of the amplitude and phase characteristic can be predicted as shown 
in figures 5(a) and (b). 

At the lower frequency range, where f < 1/(27 ttj^^), the condensing pressure re- 
sponse to input flow rate diminishes with increasing frequency in the form of a lag-lead 

with corner frequencies at l/(27TT2g) and 1/(27 tt^). 

At higher frequencies, above 1/(2 ittj^^), the fluctua- 
tions of the flow rate are absorbed by the gas phase 
mass capacity so that the condensing pressure varia- 
tions continue to diminish as frequency increases. 

It should be emphasized that this discussion ap- 
plies to the condenser with a constant interface 
(L^ = 0). When this is not the case, such as with a re- 
ceiver attached to the exit of the tube, the condensing 
pressure response to flow rate variation may change 
significantly. 

To understand how the interface may behave, the 

second of the two derived models, as shown by the 

block diagram of figure 4, will be analyzed. When it 

Figure 5. - Shape of frequency response of js assumed that all the input disturbances in equation 
condensing pressure to vapor flow assum- 
ing no interface movement (Lc • 0). (24) except the condensing pressure variations are 
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zero, the interface movement becomes related to the condensing pressure by the transfer 
function of equation (2 5) . 

In general, this represents a second-order dynamic response. In this case, however, 
the parameters are not as simply analyzed as for the previous transfer function. No spe- 
cific dynamic characteristic, that is, underdamped or overdamped, can be ascribed to 
this relation. Some trends in the response may be noted, however, by examining the 
standard natural frequency f^^ and damping ratio u of the second- order transfer func- 
tion. These are expressed in terms of a and b 


fn- 


27 t yfsi 


(32) 


a = 


b 

2 Va 


(33) 


Typical values of the natural frequency and the damping ratio for the tube used in the ex- 
perimantal study are f^ = 0. 422 and a = 2. 36. Hence, the interface by itself responds 
as a highly over damped second- order system. 

Equations (32) and (33) can be written in more basic terms as 


2tt 



cr 


3q, 


2 9w, 



Pl 

LOq/9Lc) 


(32a) 


(33a) 


From the last form several relations between the natural frequency and damping and the 
basic system properties become evident. The natural frequency of the interface move- 
ment increases and the damping decreases as the vapor pressure drop gradient 9q^9L^ 
becomes larger. A larger liquid length representing more inertia, decreases the 
natural frequency as well as the damping. The damping ratio a increases with increas- 
ing liquid flow resistance 9qj^/9wj^ and larger tube area A^. 

The analysis of the responses of the separate elements of the condenser model still 
leaves the problem of their combination and interaction. The resultant complete conden- 
ser model will now be examined. 

The mathematical form for the simultaneous variation of flow rates, interface, and 
condensing pressure is arrived at by the combination of equations (12) and (24). Concep- 
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Figure 6. - Block diagram of mathematical model for dynamic be- 
havior of condenser. 


with the dimensionless gain variables 


tually, this is the same as combining the 
block diagrams of figures 3 and 4 in the form 
of a multiple- loop feedback system as shown 
in figure 6. 

The condensing pressure is solved for 
by substituting equation (24) for in equa- 
tion (12). The result is given by 

Pc = KlpGlp(®)WH ^ K2pG2p(s)C„ 

K3pG3p(s)P^ (34) 


K 


lp‘ 


Kt(l + /3g) 


(35a) 
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(35b) 


(35c) 


The normalized transfer functions are 

Glp = + NjS + l) 

^2p ~ G^(as^ + bs + 1) 


(35d) 

(35e) 
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The definitions of the more basic dimensionless gains and time constants are 

Ta ^Pc 
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Similarly, when equation (12) is substituted into equation (24) and solved for the in- 
terface movement, 

= KjlGjlWWh + K2i,G2i,(s)C„ - K3j,G3l(s)P^ (36) 

The dimensionless gain constants are 


^1L = 


" KjKvd + ^g) 


^2L = 


KtK^a V 


(37a) 


(37b) 
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The normalized transfer functions are 


where 
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Equations (34) and (36) represent the desired mathematical model of the condensing 
process. The first equation relates the percent variation of condensing pressure to per- 
cent changes of the input variables consisting of the input flow, the coolant temperature, 
and the receiver pressure. The second equation relates the percent variation of interface 
for the same input variables. The normalized forms which are shown separate the 
steady-state gains (the K’s) from the dynamic elements (the G^s). The K's define the 
percentage change of the output variables for percentage change in the input variables in 
the steady state. The G*s are equal to unity for steady state, but define the transient or 
frequency response behavior of the condensing pressure and interface when any of the in- 
put variables change with time. 
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Figure 7. - Experimental rig for mercury condenser dynamic studies. 






The complexity of the algebra does not allow a general discussion of the dynamic 
properties of the transfer functions. However, characteristic patterns in the dynamic 
response will be noted in the next section where the frequency responses for the experi- 
mental condenser tube are computed. 


COMPARISON OF THEORY WITH EXPERIMENTAL DATA 

The many simplifying assumptions and linearizations that are made in deriving the 
mathematical model of the condensing process raise the question of how well the model 
predicts the behavior of the real system. An estimate of the accuracy of the model can 
be obtained by comparing the analytic and the experimental frequency responses. The 
tube was equipped to measure the frequency response of only two of the derived transfer 
functions (G^p and so that a comparison for only these can be made. The other 

transfer functions are not analyzed. However, their frequency response can be easily 
computed by using the constants of tables H and IV (pp- 46 and 50) . 

The first part of this section describes how the data were obtained and how the corre- 
sponding theoretical frequency responses are calculated. The results are presented in 
graphical form for two different operating points. The areas of agreement and disagree- 
ment, as well as the possible sources of error, are discussed. In the last part, an at- 
tempt is made to estimate the effect of inaccuracies of various parameters on the pre- 
dictability of the frequency responses. 


Procedure for Obtaining Frequency Response 

The test facility is shown schematically in figure 7. It consists of a horizontal 
stainless-steel condensing tube that is 90 inches long and has a 3/8-inch nominal outside 
diameter and a 1/32-inch wall. The mercury vapor, near a quality of one or slightly 
superheated, is generated by an electrically heated tube-type boiler. A large plenum is 
provided at the exit of the boiler to assure constant pressure at the inlet of the choked 
flow control valve during the induced flow variations. The condensing tube is cooled by 
nitrogen gas flowing across the tube. The gas emerges from closely spaced holes in the 
manifolds located above and below the tube. The condensed liquid exhausts through the 
exit valve into the controlled-pressure receiver. 

The tube is instrumented to record the dynamic variations of condensing inlet pres- 
sure, receiver pressure, liquid pressure drop across the exit valve, interface move- 
ment, and mercury vapor flow rate. The pressures are measured by closely coupled 
fast response electromechanical pressure transducers. The interface movement is mea- 
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sured by special means based on the principle of the mercury switch. Conducting pins, 
electrically insulated from the tube wall, are placed at 1-inch intervals along the axis at 
the top of the tube. Those pins that are in contact with the liquid mercury send **on” 
signals which are summed to produce a voltage proportional to the length of the liquid 
column inside the tube. The mercury vapor flow is determined from measurements of 
the AP across a calibrated venturi at the inlet of the tube. 

The experimental procedure in obtaining the frequency response data is as follows: 
The mercury flow rate to the condenser is set at some desired operating point by control 
of the inlet vapor flow control valve. The nitrogen gas coolant rate is kept constant by 
setting a desired manifold pressure. The liquid-vapor interface is then adjusted to the 
desired length by control of the receiver pressure. The steady-state gain is established 
by measuring the pressure and interface positions at the upper and lower limits of the 
flow variations to be used. These limits are set about ±10 percent of the steady-state 
flow Wjj. The instrumentation and recording equipment is checked and calibrated if nec- 
essary at this time. After this standard preliminary procedure, the frequency response 
tests are conducted. The vapor flow control valve is moved sinusoidally at various fre- 
quencies to introduce sinusoidal variations in mercury vapor flowing into the condensing 
tube. (It should be noted here that the valve is designed to provide linear variations of 
flow area against valve position so that in combination with the choked condition and con- 
stant upstream pressure, a sinusoidal valve movement produces corresponding vapor 
flow variations.) The resultant changes in the primary dependent variables - the con- 
densing inlet pressure and the vapor liquid interface - together with the flow variations 
are recorded simultaneously on a direct-inking recorder. The resultant sinusoidal 
traces are then used to calculate the amplitudes ratio and phase shifts at each frequency. 
The amplitude is normalized by the amplitude obtained during steady- state calibration. 

These tests are carried out for a number of operating points, that is, various levels 
of coolant flow rates, receiver pressure, and vapor flow rate. Only two examples, 
which will be called cases I and n, are presented to illustrate two distinctly different 
frequency responses. Case I corresponds to a long condensing length and a large pres- 
sure drop in the vapor region. Case n represents a short condensing length and a small 
pressure drop between inlet and receiver pressure. 

Before presenting the results, however, the method of calculating the theoretical 
frequency response should be discussed. As mentioned previously, only the two transfer 
functions and as given by equations (35d) and (37d), are considered in the 

comparison. The parameters entering into these transfer functions are summarized in 
tables n, in, and IV of appendix B. An example of the calculation of these parameters 
is also presented therein. 

Some of the required properties of the process, in particular the pressure drop in 
the liquid leg of the condenser, are not accurately known. Hence, the estimated mean 
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values are used for the calculation. The effect of Inaccuracies and changes of this and 
other uncertain variables on the predicted frequency response will be discussed in more 
detail in the Sensitivity of Predicted Frequency Response to Parameter Changes section. 

Substituting the numerical values of the parameters as given in table II into equations 
(35d) and (37 d) gives the theoretical expression for G^p and used in the compari- 

son with the experimental results. The results for case I are 


^c _ Q (g) ^ (1 + 0. 851 s)(l + 1. 48 s + 0. 075 s^) 

1 + 3. 86 s + 3. 40 s^ + 0. 208 s^ + 0. 00203 s'^ 


(38) 


c _ 


Kil^h 


= Gil(s) = 


1 + 1. 66 s + 0. 0081 s^ - 0. 000035 s^ 


1 + 3. 86 s + 3. 40 s^ + 0. 208 s^ + 0. 00203 s'^ 


(39) 


The form used for plotting the frequency responses is obtained by replacing s by the 
complex operator 27rjf. The plots and their comparisons are presented next. 


Results and Comparison 

Figure 8(a) shows the normalized amplitude and phase of the condensing pressure 
response G^p as a function of the frequency of the vapor flow disturbances for case I. 
The experimental points are shown by squares, while the theory is represented by the 
continuous line. 

The agreement in amplitude is very good, both in the shape of the curve and the ac- 
tual amplitudes over the whole frequency range. The predicted phase lag also agrees 
very well up to 0. 1 cps, but becomes less than the data indicate in the frequency range 
of 0. 1 to 1 cps. The phase data above 1 cps could not be measured for this case due to 
recording of the sinusoids with the paper speed too slow. Several factors may be re- 
sponsible for this discrepancy in the 0. 1 to 1 cps range. There may be some delay in 
the response of pressure due to the time required for establishment of thermal equilib- 
rium between the liquid and gas phases, which has been neglected in the formulation of 
the Clausius- Clapeyron equation. Another delay may be introduced by the transport time 
required for the droplets to reach the interface. This effect will be discussed later. 
There is also a possibility of error in the reduction of the experimental data due to the 
somewhat distorted waveform of the recorded pressure oscillations. 

The basic dynamic relation for case I is a lag- lead response at the lower frequencies 
corresponding to the time constants of the wall energy storage and the liquid inventory 
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(b) Case II. Short condensing lengths. 

Figure 8. - Comparison of experimental and theoretical condensing inlet pres- 
sure frequency response to vapor flow. 

capacity. At the higher frequencies above 1 cps the vapor volume capacitance begins to 
absorb the vapor flow oscillations, thus further diminishing the pressure response. 

The condensing pressure does not necessarily respond in the same manner under all 
operating conditions. This condition is shown by case n. Besides the differences in con- 
densing length noted previously, the condensing pressure and the coefficients of heat 
transfer for both the condensing and coolant sides are considerably higher for this case. 
Consequently, the parameters entering into the transfer function as shown in table II 
differ significantly for the two cases. Both the theoretical and experimental frequency 
response show a corresponding difference as can be seen in figure 8(b). 

Ag ain relatively good agreement exists between the theoretical and experimental re- 
sults. At low frequency, the pressure oscillations are small and increase rapidly with 
increasing frequency until a maximum is reached at about 0. 1 cps. The theory predicts 
an increase in amplitude by a factor of 4, while the data show an increase of about 2. 5. 
The shape of both the amplitude and phase agree well over the whole frequency range. 

The change in the pressure variations for the two operating points may be under- 
stood more clearly by comparing the relative amplitudes of the frequency responses, as 
shown in figure 9. The amplitude for case II is normalized relative to the steady- state 
variations obtained for case I. 
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Condenser 



Figure 9. •- Change in condensing pressure frequency response illustrating 
effect of interface movement. 


The difference in amplitude may be attributed primarily to the change of the sensi- 
tivity of condensing pressure to interface movement. First the percent change in heat- 
transfer area per inch of condensing length is greater for the short condensing length of 
case II than for the long condensing length of case I. Also, for case n, the marked de- 
crease in pressure drop for the shorter condensing length allows the condensing pressure 
to deviate less from the constant receiver pressure. 

As a further ramification of the effect of interface movement on the pressure re- 
sponse, the hypothetical predicted pressure response with no interface movement is 
shown as the dotted line. This curve was obtained by substituting the appropriate param- 
eter values of case I into equation (26). Note that it represents a specific instance of the 
general shape of the predicted frequency response of figure 5(a) (p. 17), which was dis- 
cussed in the theory section. 

With no interface movement, the pressure amplitudes as seen in the figure are three 
times higher at the lower frequencies than for case I. At the higher frequencies the dif- 
ferences between all three curves becomes much smaller, since the interface response 
declines above 0. 1 cps, as will be seen next. 

The frequency responses for the interface Gjl both cases are shown in fig- 
ure 10. The data points correspond to case I. The data for case n are not available due 
to lack of instrumentation at the short condensing length. 

The predicted and measured amplitudes of the interface movements agree remark- 
ably well. The phase does not agree, however, and the measured phase lag increases 
much more rapidly than the theory would predict. This discrepancy is most probably due 
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Frequency, f, cps 


Figure 10. - Experimental and theoretical frequency response of interface move- 
ment to vapor fiow for cases I and II. 

to neglecting the transport delay of the condensate liquid droplets reaching the interface. 
The theory, in fact, assumes that as soon as the vapor condenses, the droplets arrive at 
the interface. 

“TS 

As a matter of interest, the addition of a pure delay, of the form e , to the fre- 
quency response introduces desired phase shift without changing the amplitude. A trans- 
port delay time of about 2 seconds (t = 2) results in a phase shift characteristic as shown 
by the dotted line in figure 10, which agrees very closely with the measured values over 
all frequencies. 

Unfortunately there does not seem to be a ready theory for predicting this particular 
delay, since the droplets form at various positions along the whole length of the tube wall 
and it is difficult if not impossible to predict their velocity distribution. 

The dynamics of the interface can be characterized by a second-order lag and a 
first-order lead term. For a damping ratio of greater than 1 the lead term cancels with 
one of the lags so that a basic first-order lag response results as shown by case I. 

Where the system is slightly underdamped, a slight lead effect may result as shown by 
case II. The general characteristic remains that of a first-order lag. An additional lag 
comes into effect at the higher frequencies, but this is not felt until the amplitude has at- 
tenuated by a factor of 10. The major lag time constant is a function of the wall thermal 
capacitance and the liquid fill time of the condenser. At the higher frequency, the liquid 
column inertia and resistance may come into play. As discussed previously, the phase 
does not show the appropriate characteristic of a second-order system, but seems rather 
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to include a transport delay time which adds additional phase lag proportional to fre- 
quency. 

It should be noted that on some of the subsequent data runs the interface movement 
showed considerable deviation from an idealized sinusoidal response, including some er- 
ratic behavior. This condition may be due to wetting of the tube by the mercury and/or 
to formation of large slugs of mercury in the condensing region which coalesce with the 
interface in random fashion. This type of deviation from the idealized mobile column of 
liquid cannot be predicted by the simplified model. Nevertheless, examination of the fre- 
quency responses for the previous conditions indicates that the break frequencies of both 
the interface and pressure responses can still be predicted by the theory with good accu- 
racy. 


Sensitivity of Predicted Frequency Response to Parameter Changes 

Some of the constants used in the calculations of the frequency response are gener- 
ally not accurately known and may change during the operation of the process. The prob- 
lem, which is to estimate how these inaccuracies effect the success of predicting dynamic 
behavior, is studied by varying the more uncertain properties and calculating the resul- 
tant changes in the theoretical frequency responses. 

Four parameters are investigated. Two of them are the partial derivatives of the 
vapor pressure drop functions q^ (9q^/9£^ and 9q^/9Wjj). These are difficult to pre- 
dict or measure accurately, and they may have a significant effect on the way the conden- 
ser responds. The third parameter is the partial derivative of the pressure drop in the 
liquid region (9qj^/9wj^), which enters into the determination of the liquid column move- 
ment. The fourth is the steady- state tube wall temperature T^, which is used to calcu- 
late the coefficients of heat transfer and hence may introduce uncertainty in the thermal 
time constant of the wall. 

Each of the previous four variables is varied one at a time about the reference value 
of case I. The amount of variation depends on the uncertainty of the particular property. 
For each variation the transfer function parameters are calculated and the resultant fre- 
quency response compared to that of case I. The results are summarized in table I. 

The first column identifies the variable and its units. The second column shows the high 
and low values of each of the variables. The third column indicates the calculated change 
in the amplitude and phase in frequency response of the condensing pressure. The 
amounts shown are the maximum percent change in the amplitude plot, and the maximum 
difference of phase in degrees. Where the amount of deviation is significant (greater 
than +25 percent), the actual frequency response plots are presented. Only two of the 
parameter changes have a significant effect on the pressure response amplitude. These 
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TABLE I. - COMPARISON OF CHANGES IN PREDICTED FREQUENCY RESPONSES DUE TO 


UNCERTAINTY OF VARIOUS TRANSFER FUNCTION PARAMETERS 


Parameter 

variable 

Variable values 

Maximum 
changes in 
Gip(jf) 

Variation of 
pressure 
steady- state 
gain, 

Kip 

Maximum 
changes in 
GiL(jf) 

Variation of 
pressure 
steady- state 





Amplitude, 

percent 

Phase, 

deg 

Amplitude, 

percent 

Phase, 

deg 

gain, 

Kil 

aqy 

(With 

High 

1 150 

-25 

10 

17 percent 

-10 

— 

- 6. 4 percent 


ft 

Reference 

965 

Significant 

sensitivity 

1. 35 

Insignificant 

sensitivity 

0.732 



Low 

720 

50 

-10 

-30 percent 

10 

— - 

10 percent 

9qv 

(ibAth 

High 

-70 000 

25 

-5 

-15 percent 

-25 

-5 

5. 2 percent 

3Wh 

' (Ib/sec) 

Reference 

-46 500 

Insignificant 

sensitivity 

1. 35 

Insignificant 

sensitivity 

0.732 



Low 

-20 000 

-16 

5 

18. 5 percent 

22 

10 

-6. 5 percent 


(ib/ft2) 

High 

3 660 

10 

5 

1 percent 

-30 

12 

- 1 percent 


(Ib/sec) 

Reference 

2 740 

Insignificant 

sensitivity 

1.35 

Significant 

sensitivity 

0.732 



Low 

1 830 

-17 

-5 

- 1 percent 

50 

-12 

1 percent 



High 

650 

-50 

16 

No change 

-20 

5 

No change 



Reference 

630 

Significant 

sensitivity 

L35 

Insignificant 

sensitivity 

0.732 



Low 

610 

50 

-10 

No change 

10 

-5 

No change 


are the first parameter 9q^/9L^, which is defined as the partial of vapor pressure drop 
to changes in condensing length, and the fourth parameter T^, which is the wall temper- 
ature. The resultant variations in frequency responses are given in figures 11(a) and 
(b), respectively. 

The lower value of Sq^/ 9L^ has the effect of increasing the predicted frequency 
range for which a steady amplitude is maintained. Increasing brings about a 

faster decline of amplitude with frequency. 

Changes in wall temperature show a slightly different effect. Apparently only the 
time constant of the wall is changed significantly by this variation. Since the factor 
1 + T^s enters directly into the numerator of the pressure response transfer function, 
this element acts directly as a lead. As increases, the lead effect is felt earlier 
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(c) Due to uncertainty of liquid pressure drop gradient. 


Figure 11. - Variations of predicted frequency response of condenser variables to vapor flow. 
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and the pressure amplitude remains higher. Decreasing (increasing T^) allows a 
greater decline of amplitude before the lead effect is felt. This shift in corner frequency 
of the lead element is also substantiated by the shift in the phase lag. 

The fourth column of table I gives the predicted dimensionless amplitude gain of the 
pressure response under steady-state conditions, as defined by K 2 p of equation (35a). 
Only the first and second variables show a significant effect on this amplitude. 

The next column shows the change in amplitude and phase of the interface response, 
again as referred to case I. Although there is some variation in frequency response in 
all cases, only variations in the liquid pressure drop qj^ seem to be of significance. 

The result is shown in figure 11(c). The major deviations come at the higher frequen- 
cies. This agrees with expectation since the time constants of the liquid column which 
are directly proportional to 9qj^/9wj^ are significant only at high frequencies. The last 
column, which tabulates the gains of the interface movements as defined by of 

equation (37a), shows that only the two vapor pressure drop partials are of importance in 
predicting the steady-state amplitudes. 

The general conclusions that can be drawn from the study of the sensitivity of the 
predicted frequency responses to uncertainty in the process parameters are the following: 

(1) The change of vapor pressure drop with condensing length and, in general, the 
heat-transfer coefficients of the tube have a significant effect on the predicted dynamic 
response of the condensing pressure. Thus it is important to determine these particular 
variables accurately. 

(2) The prediction of the dynamics of the interface movement is not very sensitive to 
the parameter variations. Only the liquid pressure drop characteristic as defined by 
9qj^/3wL fias a significant influence on the normalized frequency response and this only 
at the higher frequencies. 


Theory Simplification 

Although it has been shown that the dynamic behavior of the condensing pressure and 
interface can be successfully predicted, the required calculations are quite long and diffi- 
cult. Furthermore, the complexity of the parameters of the two transfer functions as 
given in table II obscures the influence of the basic time constants discussed in the ear- 
lier sections of the theory. That some of the terms in the polynomials and constants are 
insignificant is substantiated by their relatively small numerical value and by the small 
effect that variations of these parameters have on the frequency response as seen in the 
previous section. Consequently, simplification of the theory is desirable and is possible, 
as will be shown in this section. 

The method of simplification is to start with the examination of the numerical values 
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of the coefficients of G^p and as given in table n. For both cases it is evident 

that the coefficients Ng, D, and C of Gjp and and of Gjj^ are more than 

an order of magnitude smaller than the other coefficients. Whether all these are small 
for every condenser is questionable; consequently, for different configurations and con- 
densing fluids a quick check of the relative size of the coefficients should be made. For 
the condensing tube studied, however, the simplifications are apparently valid for the 
various operating points. 

In addition to neglecting some of the coefficients, the form of the remaining con- 
stants A, B, Nj, and can be further simplified by dropping those terms in their 

calculations that are small. Their relative size is determined by examining the numeri- 
cal values as given in table IV. After eliminating the terms containing 9, t^, and r^p, 
the coefficients of the transfer function polynomials can be written as 




+ 


Iw w 


1 + A 


g 


(40) 


b’ = 


iS.a - ^g) 


(41) 






(42) 


m' = 

1 - ^p 


(43) 


The resultant simplified form of the transfer functions and are 


(1 + t^s)(1+n's) 

= 

1 + A’s + B’s'^ 


(44) 
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IL 


1 + a’s + B’s^ 


(45) 
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The validity of these simplifications must be checked by comparing the resultant fre- 
quency response with the unsimplified cases. The new constants are calculated by using 
the applicable numerical values used for cases I and n in table IV. The expressions for 
the frequency responses for the simplified cases I and n are given in equations (46) to 
(49): 

Case I; 



.001 .002 . 004 . 006. 01. 02 . 04 . 06.1 .2 . 4 . 6 1 2 4 6 10 


Frequency, f, cps 
(b) Case II. 

Figure 12. - Comparison of approximation to original theory for condensing 
pressure response. 
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Phase angle, deg Amplitude (normalized), Phase angle, deg Amplitude (normalized), 




Frequency, f, cps 


(b) Case II. 

Figure 13. - Comparison of approximation to original theory for interface 
response. 


Case II: 


G'ip(w) 


\ 0.338/ \ 0.0241 / 

0. 123 \0. 123/ 


1 + 


Gil(w) = 


3f 
0. 109 


l-Hjf^:M 

0. 123 


\0. 123/ 


(48) 


(49) 


These equations are plotted in figures 12 and 13. The solid lines show the previous pre- 
diction of the unsimplified theory. As can be seen from the figures, good agreement is 
obtained with the previously predicted condensing pressure responses for frequencies up 
to approximately 3 cps. The interface amplitude agrees well up to 1 or 2 cps, but the 
phase begins to show a marked difference at even lower frequencies for both cases. 

These differences are as expected since the neglected time constants would begin to show 
their effect at these .higher frequencies. 

Thus, it can be concluded that a good estimate of the frequency response behavior of 
the condenser up to approximately 1 cps can be obtained by the simplified transfer func- 
tions. The condensing pressure response to flow rate variations is given by equation 
(44). The corresponding interface movement is given by equation (45). The coefficients 
of the new transfer functions, as defined by equations (40) to (43), are much reduced in 
complexity. Their simplified form shows that basically only the three time constants 
^w’ ^Iw’ ^r’ ^®sides the various dimensionless gains, are important in determining 
the dynamic characteristics over the specified frequency range. The meaning of these 
parameters was discussed in the earlier sections where the interface movement and pres- 
sure response were treated independently. 

The only limitations of using the simplified model are that (1) the frequency response 
is not as accurate at the higher frequencies where the neglected vapor volume capaci- 
tance and the liquid column time constant become dynamically significant, and (2) the ne- 
glected coefficients and parameters must be small. 


CONCLUDING REMARKS 

The linearized analysis of the condenser tube results in a set of transfer functions 
relating the dynamic behavior of condensing pressure and interface movement to various 
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input variables. The condensing pressure to vapor flow rate transfer function for a fixed 
interface is a lag-lead type for moderate frequencies. The time constants for this trans- 
fer function are related to the heat-transfer coefficients and the tube wall heat capacity. 

At high frequencies, an additional lag arising from the vapor volume capacitance comes 
into effect. The interface movement depends on the type of end conditions used for the 
condenser. For the constant pressure receiver, the interface responds to condensing 
pressure as a second-order system where the damping and natural frequency depends on 
the liquid column inductance, the pressure drop characteristic of both the vapor and liq- 
uid, and the liquid fill time of the condenser. 

The condensing pressure transfer function, combined with the interface transfer func- 
tion of a constant pressure receiver, was used to predict the behavior of a mercury con- 
densing tube. The predicted condensing pressure frequency response to vapor flow varia- 
tions agrees with the experimental values over a frequency range between 0 and 5 cps. 

The amplitude of the corresponding interface movement is also accurately predicted. 

The measured values of the phase of the interface, however, show a much greater lag 
than predicted. This is attributed to the delay of the liquid droplets reaching the interface 
which has been neglected in the theory. 

The condensing pressure frequency response is sensitive to changes in the operating 
point especially when these changes involve the condensing length, heat-transfer coeffi- 
cients, and some of the pressure drop characteristics. The responses are not so sensi- 
tive to many other system parameters, so that good prediction may be obtained without ‘ 
accurate knowledge of all the process properties. 

Considerable simplification in the form and length of the polynomials of the two 
transfer functions relating condensing pressure and interface position to flow rate varia- 
tions is made possible due to the insignificant magnitudes of some of the coefficients. 

The simplified model agrees very closely with the more complex theory except at the high 
frequencies above 1 cps where the inaccuracies would be expected. 

Since the theory is derived for the general case of a crossflow cooled condenser, it 
may be concluded that the model will predict dynamic response of condensing pressure 
and interface movement in other condensers of this type. Where such condensers have 
entrance and exit conditions other than treated herein (choked inlet and a constant pres- 
sure received at outlet), the derivation still may be applied by addition of the proper 
pressure-to-flow impedance functions to the basic condenser model as given by equation 
(12). Several aspects of condenser dynamics, however, still need to be analyzed. Among 
these are the dynamic interaction between tubes in condensers of multitube construction 
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and the effects of coolant dynamics, particularly in liquid- cooled counterflow or parallel- 
flow condensers. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, Februaty 14, 1966. 





APPENDIX A 


SYMBOLS 


A, B, C, D 

a’,b’ 

a, b 
^s’^s 

Cn 

Dj, Dg 


Gi 

Gs 

^IL’ ^2L’ ^3L 
°'lL 

^Ip’ ^2p> ^3p 




j 


coefficients of denominator of transfer function, see table II 
cross-sectional area of tube, 7rDjy4 

simplified coefficients of transfer function, see equations (40) and (41) 

coefficients of interface transfer function, see page 15 

coefficients of transfer function, see table IV 

coefficients of pressure transfer function, see page 10 

dimensionless coolant temperature variation, 

specific heat of tube wall 

inside and outside diameter of tube, respectively 
frequency of sinusoidal disturbance, cps 
natural frequency of second-order system, cps 
transfer function of interface movement, see equation (2 5) 
basic condensing pressure transfer function, see equation (12) 
common term in condenser transfer functions, see equation (35g) 
interface transfer functions, see equations (37d), (37 e), and (37f) 
simplified see equation (45) 

condensing pressure transfer functions, see equations (3 5d), (35e), and 
(35f) 

simplified Gjp, see equation (44) 
latent heat of vaporization 

coefficient of heat transfer between condensate and wall 
coefficient of heat transfer between coolant and wall 
mechanical equivalent of heat 
complex operator, yTi 

(ZLPc/Ta)(^T^/dp^) 
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Kv 


Kw 

( 9qy/ 3wjj)/ ( dqj^/ dWj) 

^IL’^L’^SL 

dimensionless steady- state gains of interface transfer functions, see 
equations (37a), (37b), and (37c) 

^3p 

dimensionless steady-state gains of condensing pressure transfer func 
tions, see equations (35a), (35b), and (35c) 


•'cw^’l/O'cw"’! + hnw^z) 

^2w 

1-^lw 

k 

proportionality constant for condensate liquid pressure drop, see equa 
tion (B4) 


percent variation of condensing length, A^ /S. 

V/ t-p 


condensing length 

^L 

liquid column length 


mass of mercury liquid column in condenser 

^IL’ ^2L’^3L 

coefficients of transfer function for interface movement, see equation 
(37d) 

m'il 

simplified see equation (43) 

m 

slope of condensing temperature profile 

NpNg 

coefficients of pressure response transfer function, see equation (35d) 

Ni 

simplified Nj, see equation (42) 

Pc 

percent condensing pressure variation at tube inlet, Ap^/p^ 

Pr 

percent receiver pressure variation, Ap^/p^ 

Pc 

condensing pressure 

Pf 

interface pressure 

Pr 

receiver pressure 

‘iL 

pressure drop in liquid region 

% 

pressure drop in condensing region 

R 

gas constant 

s 

Laplace operator, sec'^ 

T 

a 

average temperature difference between condensate and coolant, 

- % 
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Tc(^„) 


n 


w 

T (J ) 


t 

V, 


''L 

W„ 


W 

Wt 


H 


w. 


H 


w 

X 

e 

/'L 

Pv 


H 


condensing temperature 
condensing temperature at interface 
coolant temperature 
temperature drop in condensing region 
wall temperature 
wall temperature at interface 
time, sec 

volume of condensing region, 
volume of gas phase 
velocity of liquid column 
percent condensate flow variation, Aw^/w 
percent flow variation, Awjj/wjj 
percent liquid flow variation, AWj^/wj^ 
condensate flow rate 
fluid flow rate into condensing region 
liquid outflow rate from condenser 
quality at inlet of tube 
distance from inlet of tube 

temperature averaging coefficient, 1 - (ml /2) 
dimensionless coefficient, Kp/K^Kj. 
dimensionless coefficient, K^(l + K^)/Kp 
dimensionless coefficient, (1 + K^)/K^ 
vapor density coefficient, Pg/PyHj^J 
vapor fill time of condenser, V^^/WjjX^ 
liquid fill time of condenser, V^(pj^ - Py)/Wjj 
density of liquid condensate 
density of vapor 
density of tube wall 

time constant of liquid column movement, lj^/A^(3qj^/9Wj^) 
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time constant of heat flow into wall, + h^w^l^ 

factors of quadratic term in G^, see equations (27) and (28) 
simplified and T 2 , see equations (29) and (30) 

Tp/(1 * V 

time constant of heat flow into wall, 

damping ratio of second- order system 
increment of variation of 
value of variable at steady state 





APPENDIX B 


CALCULATION OF PARAMETERS FOR CONDENSER TRANSFER FUNCTIONS 

The two transfer functions of interest relating the condensing pressure and vapor- 
liquid interface movement to vapor flow rate variations are 


^ ^ Kjpd H. T^ S)(NaS=^tNjS^l) 

Ds^ + Cs^ + Bs^ + As + 1 

^ ^ Kjl (l + 

Ds*^ + Cs^ + Bs^ + As + 1 


(Bl) 


(B2) 


The transfer function parameters and their definitions are tabulated for the various 
cases investigated as given in table n. These constants were calculated by using the 
values of the process variables and dimensions at the operating point. An example of 
these calculations is presented herein. 

The basic system constants needed for the calculation of the transfer function param- 
eters are given in table HI. The first group is easily obtained from the dimensions of 
the tube, the gage readings, or from the material properties. The second and third 
groups of constants are calculated from the steady- state data. The calculations for 
case I are presented as an example. 

The slope m is calculated directly from the condensing temperature profile as 
shown in figure 14: 


T^(0) - T^(t^) 
Tc«»‘c 


j $0 - 625 ^ (j QJ3 
680 X 6.23 


The parameter Zj^, as defined in equation (12b), gives 

2, 0.96 

^2 2 


The average temperature difference is then 

Ta = T^.(0)Zj^ - Tj^(O) = 680(0. 96) - 170 = 483° R 
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TABLE n. - DIMENSIONLESS TRANSFER 


TABLE ni. - BASIC SYSTEM 


FUNCTION GAINS AND TIME 
CONSTANT PARAMETERS 


Transfer function 

Case 

parameter 



I 

n 


(1 + 

1.35 

0. 264 

<‘ ■ V"S 

IL K^K^(l + ^g) 

0.732 

0. 948 

r = — — , sec 

4(h^^D, + h^^D2) 

0.851 

0.471 

,, _b + ^v<^ip + e) 

1 1 + ^ 

1. 485 

5. 71 

Nj. ,sec2 

0. 075 

0. 46 

^ ■ 1 + /3g ’ 

3.86 

1. 66 


3. 40 

1.73 

^ 3 

C = , sec 

1 + ^0. 
g 

0. 208 

0. 149 

s 4 

D - , sec 

0. 00202 

0. 00104 

V - + ^Ip) 

M.t = , sec 

IL j . 

1. 66 

1. 11 

M <S*VlA ^^2 ■ 

2L 

■0.0081 - 

-0. 00235 

,, ^s^lp'^p 3 

Mqt - , sec 

3L 

0.000035 

0. 00002 


CONSTANTS FOR TWO 
OPERATING POINTS 


Variable 

Case 


I 

n 

Directly read variables and constants 

Dj, ft 

0. 025J 

0.00255 

Dg, ft 

0.030( 

0.0306 

Ac. 

0.00050J 

0. 000508 

c^, Btu/ab)(°F) 

0. 11 

0. 11 

P„. Ib/ft® 

490. C 

490.0 

75v. Ib/ft® 

0.22E 

0. 33 

p^, lb/ft3 

792.0 

790.0 

Wjj, Ib/sec 

0.0315 

0. 035 

h- ft 

6. 23 

4. 16 

‘l. ft 

2. 27 

4. 34 

Pc> 

2232. 0 

3040. 0 

Tc(0), °F 

680.0 

710.0 

Xq, dimensionless 

1.0 

1.0 

Btu/lb 

126. 1 

126.0 

T„(0), °R 

630.0 

630.0 

Tw - T„. °R 

410.0 

495.0 

To - T«,. °R 

65.0 

45.0 

Calculated constants 


m, ft ^ 

0,013 

0.002 

Zj^, dimensionless 

0.96 

0.996 

T^, °F 

482.0 

537.0 

Btu 

^nw* 9 „ 

(fO<°R)(sec) 

0.0153 

0. 0224 

. Btu 

^cw’ 2 o 

(ft'‘)(°R)(sec) 

0. 159 

0. 294 

Vc=Ac?c> ft^ 

0. 00317 

0. 00213 

Parameter variables 


(lb/ft2)(sec) 
9Wj^ ’ lb 

2750 

3300 

^v (Ib/ft^) 
3«c ’ ft 

965 

750 

^v {lb/ft^)(sec) 
’ lb 

-46 500 

-65 000 

'*'^0 °R 

0. 0445 

0. 035 

^Pc ’ (Ib/ft^) 

Ts, °R 

55.0 

6.0 

T„(0). °R 

1090 

1125 

T^(0) - (Tg/2), °R 

602 

662 
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o 
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Condensing 
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80 
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100 


figure 14. - Typical temperature profile variation due to 
changes in mercury flow rate and condensing length. 


The heat-transfer coefficients h^^ and h^^ are obtained by employing the measured 
average temperature difference between the wall and coolant and wall and condensing 
vapor under steady- state conditions when the heat flux rate is known: 




nw 


- tJ 


(126. 1)(0.0315)(1.0) _ 0.0153Btu 

(3. 14) (0. 0306) (6. 23) (602 - 170) ^^^2^ 


where 


T = T (0) - 
w w' ’ 


T^(0)-T(£J cc 

5 — = 630 - — = 602 


and similarly, 


cw 




i-DjJ JTc(O) - T^(0)] 


= 0. 159 


The condensing volume is the product of the cross-sectional area and the condensing 
length: 


V = A £„ = 3. 17x10"^ ft^ 
c c c 
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Fluid flow rate, 
Wh 

2.85x10'2 

3. lxlO'2 

3.5x10"^ 


Pressure drop data 
from mapping 



Figure 15. - Vapor pressure curve 
for mercury illustrating calcula- 
tion of the local slope of the sat- 
uration curve (dT^/dpc). 


Figure 16. - Example of graphical estimation of pressure drop gradients 
dqy/dw^ and dq^/dL^ for variations about steady-state operating 
point. 


Calculation of Partial Derivatives 


The derivative dT /dp is by definition the slope of the mercury saturation curve, 
c c 

It can be calculated by using the Clausius- Clapyron relation as given by equation (7) or 
by graphical means. The later method is employed for case I as shown in figure 15: 


dT 
c 

^Pc 


0.044 



A fairly simple method was devised for evaluating the two partial derivatives 
Sq^/di^ and 3qy/9wjj. The method depends on obtaining a single experimental map of 
q against L for constant vapor flow in the vicinity of the operating point in addition to 
the three points of pressure drops against condensing length at the midpoint and the two 
extremes of the flow perturbations. The partials can then be evaluated as shown in fig- 
ure 16. The partial 9q^/8Lp is obtained from the straight line approximation to the ex- 
perimental map as indicated by circular data points. When lines of the same slope are 
drawn through the two points (indicated by a square) obtained from the flow variations and 
the vertical distance between them is measured, the other partial dq^/ 9Wjj is deter- 
mined. 
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The negative value of 9q^/ 9Wjj that results for both cases is not easily explained by 
steady-state analytic calculations. The result is substantiated, however, by the experi- 
mentally determined steady-state mapping of the condenser operation. A typical plot of 
condenser inlet pressure and receiver pressure against condensing length for a gas cool- 
ant rate, such as used for case II, is shown in figure 17. The solid lines are the con- 
densing inlet pressures with constant mercury vapor flow as the parameter. The dotted 
lines are the receiver pressures corresponding to the inlet pressures. The pressure 
difference between each pair of lines represents the pressure drop of the condenser which, 
except for the small liquid pressure drop, is equal to q^. The variation of q^ with flow 
rate can be obtained by measuring the pressure drop for each flow rate at a constant con- 
densing length. 

The last partial derivative of liquid pressure drop to flow variations 9qj^/9Wj^ is 
not easily evaluated, since the liquid pressure drop is small and the pressure transducers 
in the liquid mercury exhibit much noise. Estimation from the recorded traces indicate 
approximately 0. 2-psi pressure drop through the wide open exit valve at the midrange 
flow. Other pressure drop measurements indicate an additional liquid friction drop in 
the tube from 0. 1 to 0. 4 psi depending on the liquid length. When an orifice type restric- 
tion is assumed, the liquid pressure drop would obey a square law relation as shown in 
equation (B3) : 
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